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1. Introduction 

We continue our LO analysis 1 of the par- 
ton distributions in the unpolarized real photon, 
which led us to a construction of the CJKL par- 
ton parametrization, improving and broadening 
our investigation. Like previously our aim is to 
develop a proper description of the heavy charm- 
and bottom-quark contributions to the photon 
structure function, F^ (x,Q 2 ). The new mod- 
els FFNScj/f 1 and CJK 1 are slightly modified 
versions of our previous FFNSc jkl and CJKL 
models, respectively. In addition we analyse the 
FFNScjk 2 model, which includes the so-called 
"resolved-photon" contribution of heavy quarks 
to F^XjQ 2 ) , given by the 7*G — > hh process, 
and CJK 2 model, with an improved positivity 
constraint. All models are based on the idea of 
radiatively generated parton distributions intro- 
duced by the GRV group (see e.g. [3]). 

In this analysis we neglect TPC27 data, as in 
other recent analyses, and slightly modify, with 
respect to previous analysis, both types of mod- 
els. Moreover, for the very first time for the 
photon, we estimate uncertainties of the parton 
distributions due to the experimental errors of 
data. Following the analyses of this type per- 



formed for proton structure by the CTEQ Col- 
laboration, [U-El and tne MRST group, [7j we 
use the Hessian method to obtain sets of parton 
parametrizations allowing, along with the parton 
parametrization of the best fit, to calculate the 
best estimate and uncertainty of any observable 
depending on the parton densities in real photon. 
We compare our results for FT? c and prompt pho- 
ton production in 77 obtained with the Hessian 
and Lagrange (@],[7|, |Ej) methods. 

2. FFNSc jk and CJK models 

The difference between the FFNSc jk and CJK 
models lays in the approach to the calculation 
of the heavy quark h (c and b) contributions to 
the photon structure function F^ix^Q 2 ). The 
FFNScjif models base on a widely adopted Fixed 
Flavour Number Scheme in which there are no 
heavy quarks among partons. They contribute to 
F 2 7 0, Q 2 ) by a 'direct' (Bethe-Heitler) 7*7 -» hh 
process (F^ h {x, Q 2 )\dir), in addition one can also 
include the so-called 'resolved'-photon contribu- 
tion from 7 *G -> hh (F% >h (x, Q 2 )\ res ) M 

FZ(x, Q 2 )\ FFNS = x J2 4i.il + Q 2 ) (!) 

i=l 



2 



2 

h(=c,b) 

with (<z7) being the light u,d,s quark (anti- 
quark) densities governed by the DGLAP evolu- 
tion equations. 

The CJK models adopt the ACOT x scheme 
[2] which is a recent realization of the Variable 
Flavour Number Scheme. In this approach heavy 
quarks are, similarly to the light quarks, con- 
stituents of the photon. Therefore, apart from 
the direct and resolved-photon contributions also 
qh contribute to Fc[[x, Q 2 ) . It leads to a dou- 
ble counting of the heavy quark contributions to 
F^(x, Q 2 ) , therefore inclusion of the correspond- 
ing subtraction terms, i^Jdir.sub and i^,J rcSiSub , 
is necessary. 

Next, following the ACOT x approach, we use 
the Xh parameters to obtain the proper vanishing 
of the heavy-quark densities below the kinematic 
thresholds, for DIS given by the W 2 = Q 2 (l - 
x)/x = 4m 2 . Substitution of x with \h = x(l + 
4rn 2 t /Q 2 ) in qh and in subtraction terms should 
lead to the correct threshold behaviour as Xh ~ ¥ 1 
for W — » 2mh- In the CJK models 

3 

Fn^Q 2 )\cjK =Y / xe 2 i (q]+q7)(x,Q 2 ) (2) 

z=l 

2 

+ E xe l(Qh+€)(Xh,Q 2 ) 

h(=c,b) 

2 

+ J2 [Fl h (x,Q 2 )W + Fl h {x,Q 2 )\ res 

h{=c,b) 

2 

^ ] [-^/ildir.sub + -P^ftlres.sub • 
h(=c,b) 

The QCD evolution (DGLAP) starts from a 
scale chosen to be small, Qq — 0.25 GeV 2 , for 
both types of models, hence our parton densities 
are radiatively generated. As it is well known the 
point-like solutions of the evolution equations are 
calculable without further assumptions, while the 
hadronic parts need the input distributions. For 
this purpose we utilize the Vector Meson Domi- 



nance (VMD) model, assuming 

P had {^ Ql) = E Qo) = Ql\ (3) 

V JV J p 

where the sum over all light vector mesons (V) 
is proportional to the /? -meson parton density 
with a parameter k. We take the input densities 
of the p meson at Qq in the form of valence-like 
distributions both for the (light) quarks (v p ) and 
gluons (G p ): 

xv p (x,Q 2 Q ) = N v x a {l-xf, (4) 
xG p (x,Q 2 Q ) = N g xv p {x,Q 2 ) = N g x a {l - xf , 
xQ p {x,Ql) = 0, 

where N g = N g N v . All sea-quark distributions 
(denoted as £ p ), including s-quarks, are neglected 
at the input scale. The valence-quark and gluon 
densities satisfy the constraint representing the 
energy-momentum sum rule for p: 

[ x{2v p (x,Q 2 ) + G p {x,Q 2 ))dx = 1. (5) 
Jo 

One can impose an extra constraint related to the 
number of valence quarks: 

n v = [ 2v p (x,Q 2 )=2. (6) 
Jo 

Use of both constraints allows to express N v and 
N g in terms of a, (3 and k, reducing the number 
of free parameters to three. 

3. New analysis 

In our new analysis we improve treatment of 
the running of a s by differentiating the number 
of active quarks in the running of a s and in the 
evolution equations, and by using lower values of 
Aqc d ■ We first describe new aspects of our anal- 
ysis which are common to all considered models. 

3.1. Data 

New fits were performed using all, F%{x, Q 2 ) 
data, except the old TPC7. In our former global 
analysis £Q we used 208 F%{x, Q 2 ) experimental 
points. Now we decided to exclude the TPC27 
data (as for instance in ^H]) since these data 
are considered to be not in agreement with other 
measurements. After the exclusion we have 182 
data points. 
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3.2. a s running and values of A^ N ^ 

The running of the strong coupling constant at 
lowest order is given by the well-known formula: 



(Q 2 ) 



47T 



/? ln(Q 2 /A(^) 2 ) 



, A) 



n--N q ,(7) 



where A^' is the Aqcd value for N q quarks. 
N q increases by one whenever Q 2 reaches a heavy 
quark threshold, i.e. when Q 2 — to 2 . The con- 



ai, Nq+1 \m 2 i ) is imposed in 



dition ai^'^TO 2 ) 
order to ensure the continuity of the strong cou- 
pling constant. In the new analyses we introduce 
the number of active quarks in the photon, de- 
noted by Nf, which differs from the number of 
quarks contributing to the running of a s . 

The distinction between both numbers of 
quarks forces to use slightly more complicated 
formulae for the evolution of the parton densi- 
ties than in our previous analysis. More pre- 
cisely, we must proceed in three steps to per- 
form the DGLAP evolution. In the first step, 
describing the evolution from the input scale Qo 
to the charm-quark mass m c , the hadronic in- 
put <3 , had( x ' Qo) * s taken from the VMD model. 
In the second step we evolve the parton distribu- 
tions from m c to the beauty-quark mass, nib, a 
new hadronic input is given by the sum of already 
evolved hadronic and point-like contributions to 
the parton density. The point-like distribution 
at Q 2 = to 2 becomes zero again. The same is 
repeated for Q 2 > m 2 ,. That way we can solve 
the equations for three ranges of Q 2 , in which 
Nf = 3,4 and 5, separately. In each range values 
of bo and A depend on N q . 

In the previous work we assumed (following the 
GRV group approach [3]) that the LO and NLO A 
values for four active flavors are equal, and used 
A< 4 ) = 280 MeV from the Particle Data Group 
(PDG) report We now abandon this as- 

sumption and take for A' 4 ) = 115 MeV, which 
is obtained in LO from the world average value 
a s (M z ) = 0.117, with M z = 91.188 GeV, using 
Eq. Imposing the continuity condition for the 
strong coupling constant and m c = 1.3 GeV and 
m b = 4.3 GeV, we obtain A< 3 ) = 138 MeV and 
A® = 84 MeV. 



3.3. VDM 

In our new analysis we try to relax the con- 
straint on n v ©■ This leads to 4-parameter fit. 

3.4. FFNS 

In the FFNScjif 2 model we include the so- 
called "resolved-photon" contribution of heavy 
quarks to F^x, Q 2 ) , given by the j*G — ► hh 
process. 

3.5. Subtraction terms in CJK models 

In yQ we derived the subtraction term for a di- 
rect contribution, F^ h \dir,subtr, from the integra- 
tion up to Q 2 of a part of the DGLAP evolution 
equations, namely: 



dql(x,Q 2 ) 
dlnQ 2 



2tt 



e 2 h k(x), 



(8) 



where k(x) is the lowest order photon-quark split- 
ting function (see Eq. (7) in Ref. £Q). For the 
lower limit we took in Ref. IT] the natural for 



a Bethe-Heitler process limit Qf 



How- 



ever, since the threshold condition is W > iml, 
even for Q 2 < the heavy-quark contribu- 
tions do not vanish as long as the condition 
x < Q 2 /{Q 2 + 4to 2 ) is fulfilled. In this paper 
we take Qf ow = Q 2 , which improves quality of 
the fits. The direct subtraction term has now the 
form: 

F^ h \ dir , suH r(x,Q 2 ) = x\n%3e 4 h -(x 2 + (l~x) 2 ).(9) 

Qo t 

We apply the same change to the subtraction 
term for the resolved-photon contribution: 



p7 I 

2 h \res,subtr 



Q 2 2 a s (Q 2 ) 



Jx v y 



).(10) 



As we noticed in Section 2 the x — > Xh substitu- 
tion leads to the proper threshold behavior of all 
the heavy-quark contributions to the K7(x, Q 2 ), 
except for both subtraction terms. This is al- 
ready seen in Eq. JjJJ that this term does not 
vanish for \h ~ * 1 and therefore by subtracting 
it the resulting heavy-quark contribution to F^ 
may become negative for large x. An extra con- 
straint to avoid this unphysical situation is, thus, 
needed. In Ref. pQ we imposed the simple condi- 
tion (positivity constraint) that the heavy quark 
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contribution to F^(x, Q 2 ) has to be positive. Un- 
fortunately, this constraint was not strong enough 
and for some small windows at small and large 
x the still unphysical situation F^ h (x, Q 2 ) < 

F^ h (x, Q 2 )\direct + F2 h {x, Q 2 )\ reso ived was found 
[T2|. In this paper we apply new condition: 

F 2,h( x >Q 2 ) > Fl h {x,Q 2 )\ dlr +Fl h {x,Q 2 )\ res .(ll) 

4. Results of the new F^ global fits 

In the fits we use 182 F^{x, Q 2 ) experimental 
points with equal weights. Fits based on the least- 
squares principle (minimum of x 2 ) were done us- 
ing Minuit • Systematic and statistical errors 
on data points were added in quadrature. 

The CJK 1 and 2 models differ in the form of 
the positivity constraint. In the CJK 1 model 
we keep the old condition F2 h (x,Q 2 ) > while 
the CJK 2 model imposes the condition (JTTJ. In 
these models we do not apply the sum rule 10 and 
have four free parameters: a, /3, JV„, k related to 
the initial quark and gluon densities at the scale 
Ql = 0.25 GeV 2 ©. 

The two FFNS models differ by the resolved- 
photon contribution of heavy quarks to 
F^(x,Q 2 ). It appears only in the FFNScjk 
2 model. In both FFNS models we impose con- 
straints JHJ and JSJ. That allows to write N v in 
terms of a, f3 and k reducing the number of free 
parameters to three 1 . 

The results for the total x 2 f° r 182 points and 
the x 2 P er degree of freedom for our new fits are 
presented in table The fitted values for pa- 
rameters a, /3, k and N v are also presented to- 
gether with the errors obtained from Minos pro- 
gram with the standard requirement of Ax 2 = 1. 
In the case of the FFNScjk models the N v pa- 
rameter is calculated from the constraint (0) and 
therefore we do not state its error. Note, that 
the valence number integral n v (JBJ) gives in CJK 
models 1.94 and 2.00, for CJK 1 and CJK 2, re- 
spectively. 

The x 2 P er degree of freedom obtained in our 
new fits, between 1.5 and 1.7, is better than the 

1 The test fits without the number of valence quarks con- 
straint gave n v Kb 0.5 and rj 1.4 in the FFNScj/r 1 and 
FFNScjk 2 models, respectively. Both these values are 
very far away from the expected value 2. 



old results, mostly due to much lower as 
well as due to modification of the subtraction con- 
tributions in the CJK models. The old xVdof 
for the same set of 182 data points read 1.99 in 
the FFNScjkl and 1.80 in the CJKL model 2 . 

We observe that the only real difference in 
X 2 /dof is between the first and other three fits 
and is a result of inclusion or not of the 7*G — > hh 
contribution to F^x, Q 2 )- It is obvious that tak- 
ing this process into account improves the agree- 
ment between model and data. 

In light of these results one can conclude that in 
both, so different treatment of heavy-quark con- 
tributions to the photon structure function, all 
fitted parameters are similar. This is related to 
the fact, that we use in global fits only data for 
F2, quantity dominated by the light-quark con- 
tributions. 

4.1. Comparison with the F^ data 

The FFNScjk 2 and both CJK models pre- 
dict a much steeper behaviour of the F^{x,Q 2 ) 
at small x with respect to FFNScjk 1 fit and 
GRS LO 14 and SaSID [T5] parametrizations. 
On the other hand these curves are less steep than 
the FFNScjkl and CJKL results from |Tj. In the 
region of a; > 0.1, the behaviour of the F^x, Q 2 ) 
obtained from all fits and parametrizations apart 
from the CJK 1 model is similar. The shape of the 
CJK 1 fit at high a; is a result of the F^ h {x, Q 2 ) > 
condition which allows for F^ h (x,Q 2 ) < 

F2 Q )\direct "i~ ^2 Q )\resolved- If effects 

in the lower position of the characteristic point 
at which the charm-quark contributions to the 
F^ix^Q 2 ) appear as compared to other mod- 
els predictions. Apart from that the CJK 1 
fit gives smaller F^x^Q 2 ) values around this 
charm-quark threshold. Finally at high Q 2 and 
high x the CJK 1 model produces much lower 
structure function values than all other fits and 
parametrizations predict. We found that the 
FFNScjk 1 fit gives very similar prediction to 
the GRS LO parametrization results in the whole 
range of x. 

FigureQ]presents the predictions for F^ix, Q 2 ), 



2 Note, however that in old analysis global fits were per- 
formed for 208 points. 
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model 


X 2 (182 pts) 


X /dof 


K 




P 






FFNScjk 1 


314.0 


1.754 


9 Ofi7+0.063 
z - zo ' -0.072 


u.zoo_ 032 


7Q9+ 0189 


0.358 


5.02 


FFNScw 2 


279.8 


1.563 


9 i in+ 084 

z - iiu -0.090 


n oiq+0.054 

u.oiu_ 051 


fl ooq+0.265 
U.OZO_ 223 


0.415 


4.51 


CJK 1 


273.9 


1.539 


9 14fi+ - 154 


n 91 8+0 070 
u.zio_ 063 


0.462t°;|£ 


269+ 076 


5.00 


CJK 2 


273.7 


1.537 


i.yo^_ 124 


n 9qq+ 077 


0.8981°;^ 


4D4+ 0116 


4.93 



Table 1 



The total x 2 f° r 182 data points and per the degree of freedom, and parameters of the best fits for 
FFNScjx and CJK models. Errors are obtained from Minos with A% 2 = 1. 



averaged over medium x regions, compared with 
the recent OPAL data n °t used directly in 
our analysis. Like in our previous analysis we see 
that all FFNS type predictions (including GRS 
LO and SaSID parametrizations) are similar and 
fairly well describe the experimental data. The 
CJK models alike the CJKL model gives slightly 
better agreement with the data. 

4.2. Parton densities 

The parton densities obtained in the CJK and 
FFNScjif models are all very similar. Of course 
there arc no heavy-quark distributions in FFNS 
models. In case of the CJK models the c 7 (x, Q 2 ) 
and 6 7 (a;, Q 2 ) densities vanish not at x = 1, as for 
the GRV LO and SaSID parametrizations, but 
at the kinematical threshold, which was our aim. 
We notice that our new parton densities have all 
similar shapes but slightly higher values than the 
corresponding old CJKL distributions. In case 
of the gluon density we find that all new curves 
are much steeper than the predictions of the old 
models and the GRV LO and SaSID parametriza- 
tions, see also below. 

4.3. Comparison with c 

In Fig. |3 we present our predictions for the 
-F^o in comparison with OPAL data not 
used directly in the analysis, and results of the 
GRS LO and SaSID parametrizations. 

All our models containing the resolved-photon 
contribution (FFNScjif 2 and both CJK models) 
agree better with the low x experimental point 
than other predictions. The GRS LO and SaSID 
parametrizations also include the resolved-photon 
term but in their case the gluon density increased 



less steep than our models predict, as was already 
mentioned. Their F 2 7 c lines lie below results of 
our new fits but higher than the FFNScjk 1 
curve which is given solely by the direct Bethe- 
Heitler contribution. 

The CJK models overshoot the experimental 
point at high x while other predictions agree with 
it within its uncertainty bounds. Again the F 2 7 c 
from the CJK 1 fit vanishes at lower x than in 
other models. 

5. Uncertainties of the parton distribu- 
tions 

Following the corresponding analyses for the 
proton we consider the experimental uncertain- 
ties of the CJK parton densities. This is first 
analysis of this type for the photon, details are 
described in |18j . 

5.1. The Hessian method 

The CTEQ Collaboration g]-[H], developed 
and applied an improved method of the treat- 
ment of the experimental data errors. Later the 
same formalism has been applied by the MRST 
group in [7|. The method bases on the Hessian 
formalism and as a result one obtains a set of 
parametrizations allowing for the calculation of 
an uncertainty of any physical observable depend- 
ing on the parton densities. 

We apply the method as described in Refs. 0] 
and [S]- For the sake of clearness of our proce- 
dures we will partly repeat it here keeping the 
notation introduced by the CTEQ Collaboration. 

Let us consider a global fit to the experimen- 
tal data based on the least-squares principle per- 
formed in a model, being parametrized with a 
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Figure 1. Comparison of the recent OPAL data 
[To) for the Q 2 -dependence of the averaged over 
0.1 < x < 0.6 F2 /a with the predictions of the 
CJK and FFNScjif models. 
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Figure 2. Comparison of the structure func- 
tion F2 c (x,Q 2 )/a calculated in the CJK and 
FFNScv/f models and in GRS LO HU and SaSID 
[To] parametrizations with the OPAL data |17|. 



set of {ai,i = 1,2, •■•d} parameters. Each set 
of values of these parameters constitutes a test 
parametrization S. The set of the best values of 
parameters {a°}, corresponding to the minimal 
X 2 , X01 * s denoted as S° parametrization. In the 
Hessian method one makes a basic assumption 
that the deviation of the global fit from Xo can 
be approximated in its proximity by a quadratic 
expansion in the basis of parameters {dj} 

d d 

Ax 2 = X 2 ~xl = ^E^ite- *)^- ?)' (12) 
»=i i=i 

where Hij is an element of the Hessian matrix, 
calculated as 

^=K^) O0 - (i3) 

Since the is a symmetric matrix, it has a 
complete set of k = 1, 2...d orthonormal eigenvec- 
tors (vi)k defined by 

d 

'^ j H i j{vj) k = e k (vi)k, (14) 
3=1 

d 

^2(vi)j(vi)k = 8 jk , (15) 

i=l 



with {e k } being the corresponding eigenvalues. 
Variations around the minimum can be expressed 
in terms of the basis provided by the set of eigen- 
vectors 

d 

a i ~ a i =^2s k z k {vi)k, (16) 
fc=l 

where {z k } are new parameters describing the dis- 
placement from the best fit. The {s k } are scale 
factors introduced to normalize {z k } in such a 
way that 

Ax 2 =X 2 -X^i> 2 - ( 17 ) 

k=l 

The above equation means that the surfaces of 
constant x 2 ar e spheres in the {z k } space. That 
way the {z k } coordinates create a very useful, 
normalized basis. The (vi) k = viu matrix de- 
scribes the transformation between this new basis 
{z k } and old {ai} basis. The scaling factors s k are 
equal to y/l/e k in the ideal quadratic approxima- 
tion fEty . 

The Hessian matrix can be calculated from its 
definition in Eq. (| 1 31) - Such computation meets 
many practical problems arising from the large 
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range spanned by the eigenvalues the nu- 

merical noise and non-quadratic contributions to 
X 2 ■ The solution (an iteration procedure) has 
been given by the CTEQ Collaboration 

Having calculated the eigenvectors, eigenval- 
ues and scaling factors we can create a basis of 
the parametrizations of the parton densities, {S^ , 
k = 1, • • • , d}. Each parametrization corresponds 
to a fc-set of cii parameters defined by displace- 
ments from a° of a magnitude t "up" or "down" 
along the corresponding eigenvector direction 

a i (S±) = a° i ±t(v i ) k s k . (18) 

For each parametrization A\ 2 = t 2 . 

The best value of a physical observable X de- 
pending on the photon parton distributions is 
given as X(S a ). The uncertainty of X, for a dis- 
placement from the parton densities minimum by 
A\ 2 — T 2 (T - the tolerance parameter), can be 
calculated from a very simple expression (a mas- 
ter equation) 

a *=!(d*(^ + )-*(^)] 2 ) • as) 

Note that having calculated AX for one value of 
the tolerance parameter T we can obtain the un- 
certainty of X for any other T by simple scaling 
of AX. This way sets of {5*^} parton densities 
give us a perfect tool for studying of the uncer- 
tainties of other physical quantities. One of such 
quantities can be the parton densities themselves. 

Finally, we can calculate the uncertainties of 
the a,i parameters of the model. According to 
Eq. i(T%|) in this case ai(S^) — ai(S^) = 2t(vi)kSk 
and the master equation gives a simple expression 

Aai = T\J2v ik Sk\ . (20) 

In practice we observe the considerable devia- 
tions from the ideal quadratic approximation of 
equation l|17fl . To make an improvement we can 
adjust the scaling factors {si} either to obtain ex- 
actly Ax 2 — t 2 at Zi — t for each of the Sf sets 
or to get the best average agreement over some Zj 
range (for instance for Zi < t). Below we apply 
the second approach. 



5.1.1. Estimate of the tolerance parameter 
T for the photon densities 

We consider now the value of the tolerance pa- 
rameter T for the real-photon parton-densities 
corresponding to the allowed deviation of the 
global fit from the minimum, A\ 2 — T 2 , as de- 
scribed above. In case of an ideal analysis a stan- 
dard requirement is Ax 2 = 1- Of course this is 
not a case for a global fit to the i 7 ^ data coming 
from various experiments, and certainly T must 
be greater than 1. Unfortunately, no strict rules 
allowing for estimation of the tolerance parame- 
ter exist, as discussed in detail in |H] and [H]- We 
try to estimate the reliable T value in two ways, 
both applied to the CJK 2 fit only. 

First we examine the mutual compatibility of 
the experiments used in the fit. That gives the 
allowed Ax 2 greater than 22 and the tolerance 
parameter T ~ 5. As a second test we compare 
the results of our four fits presented in this paper 
and find the T values for which parton densi- 
ties predicted by the FFNS and CJK 1 models 
lie between the lines of uncertainties of the CJK 
2 model. These values are large due to the dif- 
ferences between the gluon densities given by the 
CJK and FFNS models. The T ~ 5 when only the 
CJK 1 and 2 gluon distributions are compared. 
For quark distributions T < 7 when we consider 
all models. Finally we estimate that the tolerance 
parameter T should lie in the range 5 ~ 10. 

5.1.2. Tests of quadratic approximation 

For each CJK model we obtained a set of the 
{( v i)k} and {sk} values, with i and k — 1, • • ■ ,4 
(since d = 4 in CJK models). We used the it- 
eration procedure from gj. Further we adjusted 
the scaling factors {sk} to improve the average 
quadratic approximation over the Zk < 5 range. 

Further we check if the quadratic approxima- 
tion on which the Hessian method relies is valid 
in the considered A% 2 range for the CJK 2 model 
(for the CJK 1 model results are similar). In the 
left panel of Fig. |3|we present the comparison of 
the x 2 dependence along each of four eigenvector 
directions (for the eigenvector k z% = 8ik) with the 
dependence of the ideal Ax 2 — z 2 curve. Only 
the line corresponding to the eigenvector 4 does 
not agree with the theoretical prediction. More- 



Figure 3. CJK 2 model. The left plot presents 
a comparison of the \ 2 dependence along each of 
the four eigenvector directions to the ideal A\ 2 — 
zf curve. In right plot analogous comparison for 
five random directions in {zi} space are shown. 
The ideal curve corresponds to A% 2 = z f = z 2 



over it has a different shape than other lines which 
results from the scaling adjustment procedure. In 
the right panel of Fig. [3] an analogous compar- 
ison for the five randomly chosen directions in 
the {zi} space is shown. For each of directions 



= z 2 and the ideal curve corresponds to 
Ax 2 = z 2 . In this case we observe greater devia- 
tions from the quadratic approximation. 

5.1.3. CJK 1 and CJK 2 parametrizations 

The errors calculated within the Hessian 
quadratic approximation are shown in tabled for 
the standard requirement of A\ 2 = 1. They 
should be compared with the slightly smaller er- 
rors calculated by Minuit and shown in table 

All parton distributions are further 
parametrized on a grid. The resulting programs 
can be found on the web-page |19| . 

5.1.4. Uncertainties of the CJK 2 parton 
densities 

In this section we discuss the uncertainties of 
the CJK 2 parton densities, the results obtained 
with CJK 1 model are very similar. 

In Figure 01 the up, strange and charm- 
quark and gluon densities calculated in 
FFNScj/f and CJK 1 models are compared 
with the CJK 2 predictions. We plot the 
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Figure 4. Parton densities calculated in 
FFNScj/f and CJK 1 models compared with the 
CJK 2 predictions. We plot for Q 2 = 10 the 
q 1 (CJK2) / q 1 (Other) ratios of parton densities. 
Solid lines show the CJK 2 fit uncertainties for 
A X 2 = 25. 



q 1 (CJK 2)/q 1 (Other models) ratios of the par- 
ton q 1 densities calculated in the CJK 2 model 
to their values obtained with other models for 
Q 2 = 10. Solid lines show the CJK 2 fit uncer- 
tainties for A\ 2 = 25 computed with the {S£} 
parametrizations. 

First we notice that predictions of all our mod- 
els for all considered parton distributions lie be- 
tween the lines of the CJK 2 uncertainties. This 
indicates that the choice of A% 2 = 25 agrees with 
the differences among our four models. We found, 
that the SaSID results differ very substantially 
from the CJK 2 ones (not shown). 

As expected the up-quark distribution is the 
one best constrained by the experimental data 
while the greatest uncertainties are for the gluon 
densities. In the case of u 1 the A% 2 = 25 band 
widens in the small x region. Alike in the case of 
s- and c-quark uncertainties it shrinks at high x. 
On contrary the gluon distributions are least con- 
strained at the region of x — > 1. All uncertainties 
become slightly smaller for higher Q 2 . 
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model 


K 


a 


/3 


N v 


CJK 1 


2 146+ 0120 

Z.11O_ 1Q7 


o-2i8±g:8g 


o.4a2±°:iSS 


2 69+ 041 
u.zoy_ Q Q35 


CJK 2 


i.yoi_ 103 


0-299±85gi 


0.898l°;?« 


n 4n4+ 066 
u.w±_ 054 



Table 2 

The parameters of the fits for CJK models with errors calculated in the Hessian quadratic approximation 
for the standard requirement of A\ 2 = 1 



5.2. Lagrange method for the uncertainties 
of the parton distributions 

The Hessian method allows the computation of 
the parton density uncertainties in a very simple 
and effective way. However, the Hessian method 
relies on the assumption of the quadratic approx- 
imation, which as we have shown in the former 
section, is not perfectly preserved. 

There exist another method called the La- 
grange multiplier method which allows to find ex- 
act uncertainties independently on the quadratic 
approximation (for the proton structure used in 
E] ' 13 and |H1). In this approach one makes a se- 
ries of fits on the quantity 

F(\,{a i })=X , ({a i }) + \X({a i }), (21) 

each with a different but fixed value of the La- 
grange multiplier A. As a result one obtains a set 
of points (x 2 (A), ^(A) which characterize the de- 
viation of the physical quantity X from its best 
value Xq, for a corresponding deviation of the 
structure function global fit from its minimum 
= X 2 (A) — Xo- In each of this constrained 
(by the A parameter) fits we find the best value 
of X and the optimal x 2 ■ F° r A = we return 
to the basic fit which gives {a- 1 } parameters and 
allows to calculate the best X value. The great 
advantage of this approach lies in the fact that 
we do not assume anything about the uncertain- 
ties. The large computer time consuming of the 
process of the whole series of minimalizations is 
a huge disadvantage of the Lagrange method. 

5.3. Examples of cross-section uncertain- 
ties in Hessian and Lagrange methods 

Finally we made a comparison of the uncer- 
tainties obtained for the CJK 2 model in Hessian 
and Lagrange methods for two physical quanti- 
ties. For the sake of limitation of the computer 



time we chose two very simple examples: F% c 
points measured by the OPAL Collaboration [T7| 
and the *yq — > 75 part of the cross-section for 
prompt photon production in 77. Results for the 
later case are presented in Fig. 5. 

6. Summary 

We enlarged and improved our previous anal- 
ysis pQ. We performed new global fits to the 
F^XjQ 2 ) data. Two additional models were 
analysed. New fits gave x 2 per degree of freedom, 
1.5-1.7, about 0.25 better than the old results. All 
features of the CJKL model, such as heavy-quark 
distributions, good description of the LEP data 
on the Q 2 dependence of the F% and on F^ c are 
preserved. We checked that the gluon densities 
of our models agree with the HI measurement of 
the G 7 distribution performed at Q 2 = 74 GeV 2 

mi, Fig. ei 3 

An analysis of the uncertainties of the CJK par- 
ton distributions due to the experimental errors 
based on the Hessian method was performed for 
the very first time for the photon. We constructed 
sets of test parametrizations for both CJK mod- 
els. They allow to computate uncertainties of any 
physical quantity depending on the real photon 
parton densities. 

Parametrization programs for all models can be 
obtained from the web-page |19j . 
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Figure 5. CJK 2 model. Lagrange and Hessian 
method results for the direct resolved (DR) part 
of the 77 — » 7g cross-section. The dashed lines 
represent the 10 to 30% deviation from the Hes- 
sian result. 
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Figure 6. Comparison of the gluon distribution 
obtained in the HI measurement performed at 
Q 2 = 74 GeV 2 H31 with the predictions of the 
CJK, FFNScjif models and GRV LO 0, GRS 
LO 55 and SaSID ^H] parametrizations with the 
OPAL measurement |17| . 



